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ABSTRACT 



Aims. Many aspects of the evolutionary phase in which Asymptotic Giant Branch stars (AGB stars) are in transition to become 
Planetary Nebulae (PNe) are still poorly understood. An important question is how the circumstellar envelopes of AGB stars switch 
from spherical symmetry to the axially symmetric structures frequently observed in PNe. In many cases there is clear evidence that 
the shaping of the circumstellar envelopes of PNe is linked to the formation of jets/collimated winds and their interaction with the 
remnant AGB envelope. Because of the short evolutionary time, objects in this phase are rare, but their identification provides valuable 
probes for testing evolutionary models. 

Methods. We have observed (sub)millimeter CO rotational transitions with the APEX telescope in a small sample of stars hosting 
high-velocity OH and water masers. These targets are supposed to have recently left the AGB, as indicated by the presence of winds 
traced by masers, with velocities larger than observed during that phase. We have carried out observations of several CO lines, ranging 
from J = 2-luptoJ = 7-6. 

Results. In IRAS 15452 - 5459 we detect a fast molecular outflow in the central region of the nebula and estimate a mass-loss rate 
between 1.2 x 10"'' M© yr"' (assuming optically thin emission) and 4.9 x 10""* M© yr"' (optically thick emission). We model the 
SED of this target taking advantage of our continuum measurement at 345 GHz to constrain the emission at long wavelengths. For a 
distance of 2.5 kpc, we estimate a luminosity of 8000 L© and a dust mass of ~0.01 M©. Through the flux in the [CII] line (158 yum), 
we calculate a total mass of about 12 Mq for the circumstellar envelope, but the line is likely affected by interstellar contamination. 

Key words. Stars: AGB and post-AGB, circumstellar matter, Submillimeter: stars 



1. Introduction 

When a low- or intermediate-mass star has completed its core 
helium burning, it goes through the AGB phase, characterized 
by heavy mass loss. Dust formation is very efficient in AGB 
outflows and about 1-2% of the mass of the envelope is con- 
densed into dust grains. After leaving the AGB, the mass-loss 
rate drops substantially and the star may become optically vis- 
ible, as its dusty shell disperses. Eventually, once it reaches a 
temperature of 20,000-30,000 K, the central star will start to 
ionize its circumstellar e nvelope (CSE) and a Planetary Nebula 
will form (lHabingL|l996l) . Dusty CSEs re-radiate absorbed stel- 
lar light and show clear signatures in the far-infrared (IR) 
spectrum. In particular, a "color-color" diagram introduced by 
Ivan der Veen & Habing (1988) has allowed for a classification 
of sources from the ratios of flux densities measured with the 
InfraRed Astronomical Satellite (IRAS) . 

The heavy mass-loss of AGB stars (with rates up to several 
10""* Mq yr"^) is one of the most important pathways for mass 
return from stars to the Interstellar Medium (ISM). 

After the AGB phase, the star "somehow" develops a less 
massive (10"^ Mq yr"') but faster (up to lO"-' km s"') wind that 
sweeps up the circumstellar material, possibly creating shocks 
and high-density shells. UV photons due to the increase in stellar 
temperature can excite, ionize and dissociate molecules, heat the 
dust, ionize the gas (when Teff > 20-30 x 10^ K). 

In the CSEs of many O-rich AGB stars, maser emission is 
detected from SiO, H2O and OH molecules. The conditions for 



the formation of maser emission from these different molecules 
are typically met in a progression of zones. The SiO masers 
arise from the closest regions to the central star, then H2O and 
OH at progressively farther distances. In stars with thicker en- 
velopes, longer periods and higher mass-loss rates - so-called 
OH/IR stars - the spectral profiles of H2O and OH masers are 
often double-peaked, and the velocity range covered by the OH 
(10-25 km s" ' ) is usually larger than that of the H2O feature 
dEngelsl l2002h . When the mass-loss process ends (i.e., at the 
end of the AGB) and the circumstellar envelope is progressively 
diluted, the maser lines may start to disappear. In a spherically 
symmetric CSE, first SiO lines vanish, then H2O, and finally OH 
(Lewis, 1998). Given these considerations, OH/IR stars with de- 
tached circumstellar shells, whose presence can be inferred from 
doubly or multiply peaked OH-maser spectra with a large veloc- 
ity spread, are good candidate objects for studying the post-AGB 
phase. 

The members of a small but growing sub-class of 
water-maser stars have been dubbed ''water-fountain" neb- 
ulae and show unique spectral and spatial characteristics 
(Likkel & Morris, 1988). The velocity spread of their H2O spec- 
tra is very large, up to several 100 km s"', and frequently larger 
than that of their OH spectra. High-resolution interferometric 
observations have shown that their masers a re in collimated jets 
with dynamical ages shorter than 150 yr (iBoboltz & Marvel 
2005). Water-maser post-AGB stars and particularly water foun- 
tains represent a very short evolutionary phase right after the 
end of the AGB; therefore they provide a unique opportunity for 
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studying the physical conditions in the molecular envelopes right 
when the shaping agent is starting its action. 

One such star is IRAS 15452 - 5459, which shows SiO 
H7O and OH masers. It was first selected bv lvan der Veen et alJ 
(1198 9) within a sample of sources detected with IRAS as a can- 
didate star in transition from the AGB to the PN phase. They re- 
quired that their targets fall within regions Illb, IV, and V of the 
van der Veen & Habing IRAS color-color diagram (i.e., evolu- 
tionary sequence from the AGB to the PN), and thus be bright at 
12 yum (> 2 Jy) and have a small IRAS variability index (< 50). 
IRAS 15452 - 5459 (hereafter 115452) falls within region IV of 
that diagram, where variable stars with very thick O-rich circum- 
stellar shells are expected. The authors assumed a stellar tem- 
perature of IC* K and a luminosity of 1.5 x lO'* L©, estimating a 
distance of about 2.5 kpc. 

f)udmaiier et al. ( 1995h performed IR K-band spectral obser- 
vations and found CO ro-vibrational lines in absorption, which 
made them conclude that the central star should have a tempera- 
ture in the range of 3000-5000 K. 

This lower stellar temperature actually appears more con- 
sistent with the fact that I15452's CSE hosts masers in all the 
species found in CSEs. Th e object was included in the s urvey 
of OH masers carried out bv lte Lintel Hekkert et al.l(ll991b . who 
detected a peculiar OH maser at 1612 MHz with 4 peaks. The 
features span roughly -20 to -100 km s"^ and their peaks are 
found around -83, -67, -47, and -33 km s"', with a central 
velocity around -57 km s"'. Subsequent observations with the 
Austraha Telescope Compact Array confirmed the detection of 
the c omponents at large r velocities and their association with the 
star (iDeacon et al.', 2004). 

IPeacon et al, (.2007.) detected the H2O maser line at 22 GHz. 
They found an irregular spectrum with many emission features, 
covering roughly -30 to -80 km s ' and pointing to a complex 
masing region. They carried out three runs of observations over 
one year and concluded that the spectral features are stable in 
velocity over a timescale of a few months, but the peak intensity 
of the individual features and the integrated spectral flux den- 
sity vary dramatically (~70%). They also report the detection of 
an SiO maser at 86 GHz with an unusual doubly -peaked profile, 
with features centered around -40 and -75 km s"'. The detec- 
tion of OH, H2O, and SiO masers and the near-IR CO absorption 
indicate that the source is an O-rich star that has very recently 
evolved off the AGB and has started its transition across the HR 
diagram. 



2. Observations and results 

2.1. /RAS 15452 - 5459 

We observed 115452 with the APEjfl 12 m telescope, which is 
located at an altitude of 5 100 m on Llano de Chajnantor, in Chile 
dGiisten e t al., 20 06). Observations were performed in both con- 
tinuum and spectral-line mode, to investigate its CSE dusty and 
molecular components. 

2.1 .1 . CO observations 

The CO observations were carried out o n 2009 October 16 and 
19. The receivers APEX-2 and FLASH (iRiscacher et al.L l2006t 



' This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fiir Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 



iHevminck et all l2006l) were used to observe the CO J = 3 - 2 
(345.796 GHz) and J = 4 - 3 (461.041 GHz) rotational lines, 
respectively. The Fast Fourier Transform Spectrometer (EFTS) 
backend available at APEX allows for a spectral coverage of 
about 1.8 GHz, split up in 8192 channels, which gives veloc- 
ity resolution s of 0.21 and 0.1 2 km s"' at 345 and 460 GHz, 
respectively ( Klein et alil2006h . 

Since the source is known from infrared imaging to be ex- 
tended over about 25", we pointed the telescope at three differ- 
ent positions. One pointing was performed at the central coor- 
dinates of the source (right ascension 15:49:11.38; declination 
-55:08:51.6), one was off-set by +7" in both RA and DEC (this 
is indicated in the following as the North-Eastem pointing) and 
a third pointing was similarly off-set by -7" (indicated in the 
following as the South-Western pointing); each off-set pointing 
was therefore about 10" away from the central coordinates. 

The observations were carried out using a position-switching 
mode, with the reference position located 1 ° away in azimuth. 
The data reduction consisted in averaging together the scans per- 
formed for each pointing and removing a low-order ( < 2) poly- 
nomial baseline. To enhance the signal-to-noise ratio, the data 
were smoothed to resolutions of 0.63 km s ' for the J = 3 - 2 
line and 2.4 km s ' for the J = 4 - 3 line. The raw spectra 
are stored in the T^ scale and can be converted to main-beam 

brightness temperature using T^\, - T'^j^, where T* is the an- 
tenna temperature corrected for atmospheric attenuation using 
the chopper-wheel method, Fe/f is the telescope forward effi- 
ciency (for APEX F,ff -0.95 at 345 GHz and 0.97 at 460 GHz) 
and Beff is the main-beam efficiency, e stimated to be 0.73 at 345 
GHz and 0.61 at 460 GHz dGiisten et a l.. 2006). The uncertainty 
of the absolute intensity scale is estimated to be about ±20%. 
Regular pointing checks were made on strong circumstellar CO 
sources and typically found to be consistent within ~3" . The 
FWHM of the main beam, d^b, is 17.3" at 345 GHz and 13.3" 
at 460 GHz. Table [T] summarizes the integration times, system 
temperatures and rms noise for each pointing. 

2.1.2. Continuum observations 

On 2007 May 28, we observed 115452 with the Large 
APEX Bolometer CAmera (LABOCA) at 870 fim (345 GHz). 
LABOCA is an array of 295 bolometers, operated in total-power 
mode at a temperature of ~ 285 mK. The noise-weighted mean 
point-source sensitivity of the array d etermined from on-s ky in- 
tegrations is 55 mJy s'^^ per channel (Siringo et al .1120091) . 

The absolute calibration is achieved by observations of pla- 
nets, namely Mars, Uranus, and Neptune. The overall calibration 
accuracy for LABOCA is about 10% (Sirin go et al., 2009) . 

A major issue in performing sub-mm observations is the at- 
mospheric opacity. Two independent methods are typically used 
for determining the sky opacity at APEX. The first one relies 
on the measurement of precipitable water vapor performed ev- 
ery minute by the APEX radiometer along the line of sight 
(iPardo et al.l l200lh . This approach is limited by the knowledge 
of the passband, the applicability of the atmospheric model and 
the accuracy of the radiometer The second method uses skydips, 
i.e. measurements of the power of the atmospheric emission as 
a function of the airmass. The values of sky opacity resulting 
from skydip analysis are robust, though up to -30% larger than 
those obtained by the radiometer This systematic difference is 
not well understood, therefore, to reconcile the results obtained 
with the two methods, a linear combination of radiometer and 
skydip values has been used, so as to obtain the most consistent 
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3.3 


460 
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83 
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0.7 


444 


213 


8 


798 
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Table 1. Summary of the observations of IRAS 15452 - 5459 performed on 2009-10-16 (J = 3 - 2) and 2009-10-19 (J = 4 - 3). 
The rms noise is calculated after smoothing the spectra to 0.63 and 2.4 km s-1 for the J = 3 - 2 and J = 4 - 3 line, respectively. 

^ The North-East pointing was offset by 7" in RA and 7" in DEC. 
'' The South- West pointing was offset by -7" in RA and -7" in DEC. 



calibrator fluxes at all elevations, as indicated by ISiringo et al.l 
(I2009h . 

The observations were carried out in raster-spiral mapping 
mode. In the basic spiral mode, spiral scans are performed at 
a constant angular speed and an increasing radius around the 
pointing position, producing fully sampled maps of the whole 
field of view. These spirals are the preferred observing mode for 
sources brighter than a few Jy. For fainter sources, the basic spi- 
ral pattern can be combined with a raster mapping mode (raster- 
spirals) on a grid of pointing positions resultin g in a denser 
samp ling of the maps and longer integration time (ISirineo et al.L 
I2009h . 

The data reduction was performed with the Bolometer ar- 
ray data Analysis (BoA) software, specifically developed for 
APEX bolometer data. BoA allows the observer to perform all 
of the necessary reduction steps: opacity correction, flux cali- 
bration, bad channel flagging, removal of correlated noise, de- 
spiking, data weight ing, creation of a map, and image analysis 
dSiringoetalL 120091) . 

The nebula shows a flux density at 345 GHz of 260 + 40 
mJy, where the error includes the rms noise and the absolute 
calibration uncertainty. 

2.2. coin water-maser stars 

Using the APEX telescope, we made an attempt to find CO ro- 
tational emission in six selected objects known to host high- 
velocity water masers, like 115452. The sample included four 
water-fountain nebulae toward which no CO emission has pre- 
viously been reported. All targets are listed in Table |2] The 
six sources were observed in the CO(2 - 1) transitions with 
the APEX-1 receiver, five of those were observed also in the 
CO(3 - 2) transitions with the APEX-2 receiver , and two of thern 
were also observed with the CHAMP-H array (iKasemann et al.L 
l2006l) giving spectra centered on the CO(6 - 5) and (7 - 6) 
lines. The observing dates and noise levels are given in Table |2] 
The telescope beam widths and main beam efficiencies at the 
observed frequencies are given in Table |3] As a back-end, we 
used the digital EFTS, which provided spectra at resolutions be- 
tween 0.2-0.6 km s ' and bandwidths of at least several hundred 
km s"'. 

T he stars were s electe d fro m currently-av ailable maser sur- 
veys (iDeacon et al.l (l2004l) and iDeacon et al.l ([2007)) as having 
doubly peaked maser profiles with velocities larger than typi- 
cal AGB values (10-25 km s '), IRAS colors typical of post- 
AGB stars, and being observable from the APEX site. All of 
them host H2O masers and are suppos e d to ha ve just approached 
their post-AGB evolution. Ilmai et aP (l2009h observed 9 of the 



13 water fountains currently known in the CO J = 3 - 2 line 
with the Atacama Submillimeter Telescope Experiment (ASTE). 
Their rms noise ranged from 4 to 62 mK and they detected the 
line in 2 targets. In o ur selection , we have excluded the water 
fountains observed bv llmai et al.l (2009) and inc luded the newl 
discovered water fountain IRAS 1 9 1 90 + 1 1 02 (.Day et 
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This leaves us with 4 water fountains and an additional 2 stars 
with water masers. 

The observations in the CO(2 - 1) and CO(3 - 2) lines were 
obtained in position switching (total power) mode. Many of 
the sources that we observed are located close to the Galactic 
plane, therefore their reference positions were chosen 1-3 ° away 
from the Galactic plane, to avoid emission at the OFF posi- 
tions. All the observed lines of sight cross interstellar molec- 
ular clouds, which give rise to prominent emission features, 
especially strong in the CO(2 - 1) line. For all sources other 
than OH 12.8 - 0.9, these interstellar emission features occur in 
the range of velocities where the maser emission is observed. 
This strongly hampers the identification of any potentially exist- 
ing CO emission from the circumstellar material in our objects. 
Indeed, for a majority of the objects, strong interstellar lines 
overlap in spectral regions close to the expected source veloc- 
ity. For IRAS 18103 - 1738, we find one feature that is centered 
around the stellar velocity, when we compare the CO and OH 
maser spectra (Fig. [1}, but the CO feature may be of interstel- 
lar origin. In the case of two other objects, IRAS 18043 - 21 16 
and IRAS 18327 - 0715, we found relatively weak emission 
features that appear close to the expected central velocities, 
i.e. at 38 km s"' and 87 km s"' for IRAS 18043 - 2116 and 
IRAS 18327 - 0715, respectively (Fig.llli. These emission fea- 
tures are candidates for the circumstellar emission but their in- 
terstellar origin cannot be excluded without careful imaging ob- 
servations and accurate positions. For these two sources, we ob- 
tained deep integrations in the CO(6 - 5) line (simultaneously 
with a spectrum in the CO(7 - 6) line). In these high rotational 
transitions the contamination from interstellar clouds of low ex- 
citation should be negligible and the circumstellar gas can be 
identified more easily, if excited to high J levels. In both cases 
the spectra do not show any emission in CO(6 - 5) nor in (7 - 6) 
at the noise levels given in Table |2] Our search for CO emission 
in these sources is therefore inconclusive. 

In the case of OH 12.8 - 0.9, the maser features appear in 
a velocity range that is not contaminated by strong interstellar 
CO emission. We detected weak and relatively broad emission 
features in CO(2 - 1) and (3 - 2) lines covering the same veloc- 
ity range as bracketed by the m aser peaks of H?P and OH , i.e. 
between -70 and -35 km s ' (iBoboltz & Marvell l2005h . The 
detection is marginal, but the feature clearly appears above the 
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Table 3. Telescope specifications for the observed frequencies. 



CO 


Frequency 


Beam" 


'7mb 


transition 


(GHz) 


(arcsec) 




2-1 


230.538 


27 


0.75 


3-2 


345.796 


17.3 


0.73 


4-3 


461.041 


13.3 


0.61 


6-5 


691.473 


9 


0.56 


7-6 


806.652 


7.7 


0.43 



Notes. *"' Half-power beam width. **' Antenna main beam efficiency. 



3cr noise level after smoothing the spectra to a resolution of 5- 
6 km s"\ as can be seen in Figure |2] The fact that it appears 
in both lines and the good match with the maser velocities sup- 
port the detection, although more sensitive observations are de- 
sirable. 
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3. IRAS 15452 - 5459: CO spectra 

We have obtained three spectra coiTesponding to the central-star 
coordinates as derived from observations performed with the 
Hubble Space Telescope and two positions 7" ofF-set in both 
right ascension and declination. These are shown in Figure |3] 

The spectra coiTesponding to the central pointing appear as 
a superposition of at least two components: a broad feature and 
a narrow one. The broad feature seems to show an asymmet- 
ric profile, its red side being relatively abruptly cut off". Such 
a shape is likely due to the presence of an absorption fea- 
ture around -20 km s"' clearly detected in the off'-set point- 
ings, which is probably masking the line profile. T he absorption 
was al ready noticed in the OH maser spectra by iDeacon et al.l 
( |2004|) . who also pointed out that a similar feature is present 
around -55 km s ' . These are likely spurious features due to the 
ISM emission in the off-source position. In Figure |4] we com- 
pare the CO f eatures with the OH maser spectrum published by 
iDeacon et alJ (2004) and notice that the wide CO component ex- 
tends over a much larger velocity range than the maser features. 
This broad component is not detected in the off'-set pointings, 
which indicates that the part of the outflow from which this fea- 
ture originates is located in the inner region of the CSE. This is 
confirmed by the detection of the broad feature in all pointings 
in the J = 4 - 3 line, which is a better tracer of the inner region. 

In the J = 3 - 2 ofF-set pointings, we detect only a central, 
strongly asymmetric double-hom component, whose red feature 
is much weaker than the blue one. At a close inspection, the 
spectral structure looks more like a four-peak feature, resem- 
bling an optically thin, resolved double shell. 

3.1. Kinematical components 

The data acquired North East and South West of the central 
star are characterized by a double-peak structure. The blue peak 
shows a double-peak substructure that is easily seen also in the 
south-western pointing. An analogous substructure is probably 
responsible for the shape of the red peak. This closely resembles 
what is observed in the OH masers, where two features centered 
around -57 km s"^ also show multi-peak substructures. 

By fitting Gaussians to the peaks, we derive for the NE fea- 
ture a central velocity of -57.5 km s"' and a velocity width 
(FWHM) of 27.4 km s ^ while for the SW feafiire we obtain 
-55.5 and 28.6 km s"^ respectively. 
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Fig. 1. Overlay of OH maser features on our CO spec- 
tra for targets with maser data ava ilable through on -line 
archives. The maser data are from Sevenster (12002 1 and 
they are available through VizieR (lOchsenbein et al.L l2000t 

http://vizieru-trasbg.fr/viz-bin/VizieR). 
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Fig. 3. Circles as wide as the APEX beam at 345 GHz are centered around our pointing coordinates: a white circle indicates the 
central pointing and two green circles indicate the two pointings 10" off-set. The I-band background image of 115452 was retrieved 
from the HST archive (North up and East left). Red bars indicate the spectra obtained for the J = 3 - 2 CO line and blue bars those 
for the J = 4 - 3 line. 
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Table 2. Characteristics of CO observations in water fountain sources. 



IRAS ID 



OH 



Coordinates 

aj2000 (5j2000 



Vlsr" 



Date 



CO 



(km s ) of observation transition mK 



18327-0715 024.692 + 00.235" 18:35:29.82 -07:13:09.4 

18043-2116 009.097-00.392" 18:07:21.20 -21:16:14.0 

15103-5754 320.9063-00.2929'' 15:14:18.47 -58:05:20.9 

19190+1102 046.2559-01.4764'' 19:21:25.30 +11:08:40.0 

012.816-00.895'^ 18:16:49.23 -18:15:01.8 

18103-1738 012.973 + 00.133" 18:13:20.24 -17:37:17.3 



42.0 


2010-07-06 


2- 


-1 


77 




2010-07-07 


3- 


-2 


19 




2010-07-09 


6- 


-5 


48 




2010-07-09 


7- 


-6 


165 


87.3 


2010-07-05,06 


2- 


-1 


79 




2010-07-07 


3- 


-2 


26 




2010-07-09 


6- 


-5 


91 




2010-07-09 


7- 


-6 


298 


47.0 


2010-07-06 


2- 


-1 


71 




2010-07-07 


3- 


-2 


22 


50.0 


2010-07-06 


2- 


-1 


53 




2010-07-07 


3- 


-2 


25 


55.8 


2010-07-05,06 


2 - 


-1 


25 




2010-07-07 


3- 


-2 


23 


17.0 


2010-11-16 


2- 


-1 


35 



Notes. '"' Object LSR radial velocity from Deacon et al. ('2004') and 'De acon et alj ( 120071) . '*' Noise rms level of the spectrum given per bin of 
1 km s"' and in the T^b units. <'> From .Sevenster. (.2002) ''" From SIMBAlJl 



As can be seen in Figure|4] the narrow peaks match very well 
with the velocity ranges of the two central featu res observed in 
the OH maser spectrum bv iDeacon et al.l (l2004l) . They are also 
centered around the velocity accredited to the central star from 
the analysis of the OH maser lines. In spite of this, such narrow 
features cannot be conclusively distinguished from ISM contam- 
ination without carefully mapping the area around our target. 

On the contrary, the broad component detected in the central 
pointing can be certainly attributed to 115452. We fit this com- 
ponent with a function of the form 



T„Av) = r7* 



1 - 






W-Vc\<Ve 



(1) 



and null elsewhere. As central velocity Vc we take -57.5 km s 



and find T 



peak 
mb 



0.8 K and y, = 45.0 km s" 



3.2. Mass-loss rate 

We do not have observations of '^CO lines that would allow us 
to estimate the optical depth of the '^CO lines. We can calculate 
the mass-loss rate associated to the outflow generating the broad 
component from the correlation between this parameter and the 
r,„i at the center of the line. 

For a spherically symmetric envelope, if we assume that its 
emis sion is optically thin, w e ca n estimate the mass- loss rate 
from lKnapp & MorrisI ( Il985h and lOloff^sonet al.1 d 19931) as 



M = 4.55x10 



-19 



Tmh 



Log(W/0.04-)s(J) 



5/6 



/c 



1 T,Il/6 

CO ^ exp 



{D9fl^ (2) 



where T^t is the line peak main-beam temperature in K, W is 
the ratio of the flux emitted by the target at 4.6 /^m and that 
emitted by a star at the same distance with radius 5 x 10'^ cm 
approximated by a Planck curve at 2000 K, s{J) is a correction 
factor to account for the specific J — > J - 1 transition, fco is the 
abundance of CO relative to H2, Vexp the expansion velocity of 
the outflow in km s ', D the distance to the star in pc, and 9 the 
beam size in arcsec. 



From our fitting to the component, we derive T,nb - 0.8 K 
and Vexp - 45 km s"'. To estimate the W parameter, we per- 
form synthetic photometry from the observed ISO spectrum of 
the source and calculate the IRAC 4.5 /im flux density with the 
Spitzer data reduction package SMARlQ (Higdon et al, 2004, 
obtai ning 1 1.38 Jy (the sourc e is included in the GLIMPSE sur- 
vey (IChurchwell et aI.L l2009l) . but its photometry at 4.5 //m is 
not available through the on-line GLIMPSE catalog, probably 
because of saturation). We thus obtain W = 0.75. The correction 
factor i(3) is taken as 0.43, following Groenewegen et al. ( 1999]) 
and fco = 3 x 10"*, which is typical of O-rich circumstellar en- 
velopes. Finally, the distance is set to 2.5 kpc and the beam is 
17.3". The derived mass loss rate is 1.2 x 10""^ M© yr"'. 

Under the assumptio n of an op tically thick line, the mass loss 
rate can be estimated as dGroenewegen et al., 1999) 



M= 1.4- 



T„i,V^D^0^ 



2 X lO'^/Of ,9(7) 



(3) 



which returns a value of 4.9 x lO"'* M© yr '. In both cases we 
obtain quite high values, but the calculations assume spherical 
symmetry, therefore their results are probably overestimated. 

Still distinguishing between optically thin and thick regimes, 
we can calculate the mass of molecular gas corresponding to 
each spectral component. The H2 masses were derived by inte- 
grating over the whole spectrum the following formula for the 
/-th channel 

Mf = MH^D^nNf^'fco (4) 

where the CO column density Nf^ is given bv lChoi et al.l (Il993h 



A^P = 1.1 X 10'^ 



Tmb^V 



T-il 



D{n, Tk) 1 - exp(-T32) 



(5) 



D{n, Tk) was shown to be equal to 1.5 within a factor of 2 in 
the range 10 < T^ < 200 K and 10"^ < n < 10'' cm ^ m^, is the 



^ SMART was developed by the IRS Team at Cornell University and 
is available through the Spitzer Science Center at Caltech. 
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Fig. 6. Top: The two CO lines detected are plotted on each other, after resampling the spectra at a resolution of 2.5 km s ' . Bottom: 
the ratio of the two lines is shown for points where the emission from both lines is above their respective noises. 



mass of a molecule of hydrogen, D the distance to the star in cm, 
and Q the beam in sr In Table |4] we list values for T32 « 1 and 
for T = 5 as a fiducial for the optically thick case, when indeed 
the mass depends on the real value of t. We also calculated the 
energy and momentum associated to the outflow along the line 
of sight by integrating P,- - M,|y, - vqI over the spectrum (the 
central velocity vq is -57.5 km s"') and then E - j^ . 

The values of momentum and energ y can be compared with 
those found by iBuiarrabal et al.l (1200 lb . They find that most of 
the standard proto-PNe show values of kinetic energy of their 
outflows in the range 10'*'*-10'*^ erg, with low-mass stars and hy- 
pergiants respectively below and above this range. Depending 
on the opacity of the features observed, we find in the molecu- 
lar envelope around 115452 kinetic energies in the range 10"^^- 
10"*"^ erg, which would then imply its classification as a low-mass 
star^ 

IBuiarrabal et al.l (12001 1) also find that if the origin of the out- 
flow is radiation pressure on grains, an upper limit to the scalar 
momentum is given by P < L/c (g cm s"' yr~^) X 2000 (yr). For 
our source, L/c = 3.24 x 10"''* g cm s"' yr"', therefore the upper 
limit to P is 6.48 x 10-"' g cm s"' . We obtain for the momentum 



6.96x 10^^ g cm s ' under the optically-thin assumption. Even in 



g cm s 
'gems"' 

this conservative case, the momentum is larger than what can be 
afforded by the radiation: only action of radiation on grains can- 
not explain th e outflow we observe in our target, as has already 
been found bv lBuiarrabal et al.l (12001 1) for many proto-PNe. 

4. IRAS 15452 - 5459: IR archive observations 

The star has been imaged with the Hubble Space Telescope 
(HST) at optical (V) and near-lR (1, J, H, and K) wavelengths 



M 


Mh, 


P 


E 


10-* Mo yr-' 


10-^Ma 10"^'' g cm s-i 


10*2 erg 




T« 1 






1.2 


2.1 

T = 5 


6.96 


5.80 


4.9 


80 


284 


254 



Table 4. Physical parameters derived for the outflow from the 
J = 3 - 2 line. 



bvl Sahaiet al.l (|2007|) . The V-band images do not show any ne- 
bulosity and only the central star is weakly detected. The nebula 
becomes brighter at longer wavelengths, exhibiting an hourglass 
shape (~ 25"), as described by Sahai et al.. (,2007,) . to which we 
refer for a morphological description. 

4. 1 . The ISO spectrum 

115452 was observed with the Short and Long Wavelength 
Spectrometers on-board the Inf rared Space Observatory (ISO). 
Its IR spectrum was shown by ' Garcia-Lario et al.1 (12000"), who 
reported the detection of the mid-lR features attributed to 
Polycyclic Aromatic Hydrocarbons (PAHs) at 6.2, 7.7, and 8.6 
jum, a plateau between 10 and 15 yum also attributed to PAHs, and 
the 30/vm feature due to MgS. All of these features are typical 
of C-rich environments. The presence of both PAH features and 
the masers would make this object a member of the rare class of 
mixed-chemistry sources, whose nature is not well understood. 
These may be stars that have recently turned their chemistry 
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Fig. 2. CO spectrum observed toward OH 12.8 - 0.9. In the bot- 
tom, the overlay of the smoothed spectrum and the OH maser 
spectrum is shown, as in Fig.lT] 



from O-rich into C-rich and during the transition show prop- 
erties of both CSEs, or they may host structures that retain O- 
rich material, in spite of the chemical change in the overall CSE. 
Another option is that photodissociation makes atomic C and 
O available for reaction and then both O-bearing and C-beari ng 
chemical species are found in the CSE (iLittle-Marenini [l986l) . 

The spectrum also shows several atomic forbidden emission 
Unes of [Olll], [Nlll], [MI], [CII], and [OI], which are typi- 
cal of photodissociation regions (PDRs), but the continuum is 
clearly contaminated by spurious emission that can be attributed 
to Galactic cirrus. We have retrieved the highly-processed data 
set from the ISO web archive and in Table |5] we list the lines 
detected in the spectrum, their widths and their fluxes. Indeed, 
as can be seen in Figure |7] the area where the source is located 
is rich in molecular gas and dust, and, since a background spec- 
trum is not available, it cannot be ruled out that this diffiise gas 
affects or is responsible for the spectral features observed, keep- 
ing in mind that the FWHM of the ISO beam was ~ 1 .5"-93" in 
the 2-245 fim range. 



• CO J=3-2 East 
" " • OH 16ia MHz (scaled) 




Fig. 4. Overlay of ou r CO J = 3 - 2 data and OH 1612 MHz 
(iDeacon et al.L llool lines for IRAS 15452 - 5459. The off'-sets 
from the central position in (RA, DEC) for the Noth-East and 
South-West pointings were (7, 7) and (-7, -7), respectively. 
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Line 



Peak Wavelength FWHM 



yum 



yum 



Flux 

10-1** Wcm-2 



[CII] (157.74) 
[OI] (145.52) 
[Nil] (121.90) 
[OIII] (88.36) 
[OI] (63.18) 
[OIII] (51.81) 
[Mil] (57.33) 
[Sill] (34.84) 
[SIII] (33.48) 



157.724 + 0.002 
145.54 + 0.03 
121.97 + 0.02 

88.367 + 0.008 

63.191+0.007 
51.79 + 0.02 
57.32 + 0.01 

34.793 + 0.003 
33.46 + 0.01 



0.662 ± 0.002 
0.65 ± 0.03 
0.99 ± 0.02 

0.316 + 0.007 
0.24 + 0.07 
0.45 + 0.02 
0.26 + 0.01 

0.09 + 0.007 

0.160 + 0.009 



5.42 + 0.04 

0.16 + 0.02 
2.0 + 0.1 
3.6 + 0.2 

0.79 + 0.07 
3.9 + 0.6 
1.7+0.2 

0.43 + 0.03 
2.4 + 0.4 



Table 5. Atomic lines detected in the ISO spectrum of 115452. In column one, the expected peak wavelengths are listed in //m. 



4.2. PDR lines and circumstellar radiation field 

In a study of fa r-IR PDR li nes in AGB stars, proto-PN e, and 
PNe. lFonget all (12001) and ICastro-Carrizo et all (1200 lb notice 
that the detection of PDR lines depends mostly on the temper- 
ature of the central star, and that shocks, massive outflows, and 
the interstellar radiation field do not seem to play a role. On the 
basis of their observations, they distinguish circumstellar from 
interstellar PDRs by the flux ratio of the [CII] (158 /jm) to [OI] 
(63 yum) lines. The authors in fact observe that the former is al- 
ways less intense than the latter, when both lines are clearly at- 
tributed to a circumstellar PDR. This clear attribution was not 
possible for their coldest sources, for which background spec- 
tra indicated interstellar contamination. This can bias the con- 
clusion that in circumstellar PDRs I[cii] can be expected to be 
smaller than I[oi] . 

To characterize the average radiation field around the star, we 
can introduce the parameter G, the intensity of the UV radiation 
field in Habing units tCastro-Carrizo et al.n2001i) : 



G = 



LF,. 



AttR^Go 



(6) 



where L is the stellar luminosity, Fuv the fraction of energy at 
wavelengths shorter than 6 eV (~ 2066 A) in a blackbody spec- 
trum at the same temperature as the star, Ri is the inner radius 
from the star to the PDR and Go = 1 .6 X 10"^ erg cm"^ s ' is the 
average UV flux in the diffuse ISM. 

The ratio between [CII] (158 fj.m) and [OI] (63 yum) lines 
depends on the physical conditions in the PDR and the [CII] 
line is known to be the dominating coolant for average radiation 
fields with G < 10^, while the [OI] dominates at higher values of 
G. Therefore, although contamination is certainly possible, the 
atomic lines observed in the ISO spectrum of 115452 may indeed 
arise in a circumstellar PDR with a weak average radiation field. 

We can calculate what value of G we expect for a star like 
115452. The estimation of/?, can be done by SED modeling, tak- 
ing as Ri the inner radius of the dust shell. To model the SED, 
we have collected phot ometric data from s everal catalogs and 
used the code DUSTY dlvezic et all 1 19991) to solve the equa- 
tion of radiation transfer in the circumstellar dust shell. DUSTY 
only treats spherical symmetry, which is clearly not our case, 
but it can still provide us with an approximate estimate of the 
inner radius. In Table |6] we have listed the data collected from 
on-Une catalogs, namely USNOBl, 2MASS, GLIMPSE, WISE, 
AKARI and IRAS; the data are shown in Figure [8] 

We ran our model with T* - 3500 K, standard density dis- 
tribution of r"^, grain size distribution with minimum radius 
flmm - 0.005 and maximum a„,ax - 0.25 yum varying as a""-'^ 



and optical constants for amorphous silicates (appropriate for 
O-rich environments). With the exception of the mid-IR range, 
where the excess of emission can be regarded as an indication for 
a disk, a good fit to the observational data can be obtained with a 
distance of 2.5 kpc, a luminosity of 8000 L©, and a circumstellar 
absorption Ay = 58 (r ~ 1 around 3.35 /im). 

From our DUSTY modeling, we can estimate a FIR con- 
tinuum flux of 1.7 X 10"" W m"^ and an inner radius of 
5.7x10'^ cm. Therefore we findG ~ 170,withFf/y = 3.4x10"''. 
115452 seems then to fall within the regime where the [CII] line 
is the dominant coolant. 

To estimate the intensity of the FUV field and the density 
in the PDR, we use the models described in iKaufman et alJ 
(l2006h and available on-line at http://dustem.astro.umd.edu/ For 
a given set of gas-phase elemental abundances and grain prop- 
erties, each model is described by a constant H nucleus density, 
«, and incident far-ultraviolet intensity G. The models solve for 
the equilibrium chemistry, thermal balance, and radiation trans- 
fer through a PDR layer The values of rin, and G are obtained 
by overlaying two sets of contour plots and finding the intersec- 
tion of the observations. The procedure requires either two sets 
of line-intensity ratios or one set of line-intensity ratios and the 
ratio of the [CII] + [OI] intensity to that of the infrared dust con- 
tinuum emission. The abundances used in the calculation were 
[C/H] = 1.4 X 10-4, [O/H] = 3 X 10-\ [Si/H] = 1.7 x 10"^ 
[S/H] = 2.8 X 10"^ [Fe/H] = 1.7 x 10"^ [Mg/H] = 1.1 x lO"'' 
and [PAH/H] = 4 x lOl 

Either a filling-factor correction should be applied to each 
line or line intensities should be used, as lines of different species 
do not necessarily arise from the same regions. For example, 
diffuse [CII] emission may fill the observing beam, whereas [OI] 
emission may arise from a smaller region at higher density and 
temperature. A similar correction should also be applied to the 
FIR continuum. 

In Figure |9] we show contour plots for the ratios between the 
sum of the two PDR lines and the FIR continuum as a function 
of G and n, obtained from the PDR modeling. In the rough as- 
sumption that all of the emitting components arise from regions 
with similar sizes, we have (/[o/j + Ac//])/^f/« ~ 0.0037. For our 
calculated value of G, this corresponds in Figure |9]to a density 
of 1 17 cm"^. This density value is somewhat low for a post-AGB 
CSE, but may be found in the tenuous region where the fast wind 
has swept away most of the material. 

4.3. Envelope mass 

We can derive the mass of dust by simply fitting a modified 
blackbody to the observational data beyond 20 yum, neglecting 
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USNOBl 


A [jt/m] 
F, [mJy] 


0.7 
0.042 + 0.004 


0.9 
0.90 + 0.09 






2MASS 


F, [mJy] 


1.235 
61.5 + 8.7 


1.662 
299 ± 56 


2.159 
1080 + 30 




IRAC 


A\pm] 
Pa [Jy] 


3.6 
3.470 + 0.49 


5.8 
18.52 + 0.76 






WISE 


A [jum] 
F, [Jy] 


3.35 
3.26 + 0.37 


4.6 
15.81+0.73 


11.56 
82.77 + 0.61 


22.24 
201.8 + 3.7 


AKARI 


A\pm] 
Fa [Jy] 


9.0 

44.8 + 3.9 


18.0 
171.4+10.7 






IRAS 


A\pm] 
Fa [Jy] 


12 
87 + 9 


25 
243 + 12 


60 

274 + 63 


100 
401 



Table 6. Photometric data on 115452 collected from on-line catalogs. Infrared data are available on line through the IRS A archive 
(http://irsa.i pac.caltech.edu ), while optical data are available through VizieR tOchsenbein et al., , i2000i, http://vizier.u-trasbg.fr/viz- 
bm/VizieR). 



the IRAS point at 100 //m, because clearly contaminated. This 
procedure provides us with an average dust temperature Td and 
a dust emissivity index P that can be used to calculate the mass 
of dust following lBeuther et al (|2005|) : 

2.0 X 10-4 Fy ( D \^ ( V \-3-/* 



Jv(Td) 



hv\ 



n-l 



The formula assumes an average grain size of 0. 1 //m and 
grain density of 3 g cm"^; h and k are the Planck and Boltzmann 
constants, D is the distance to the star (2.5 kpc), Fy the flux den- 
sity measured at the frequency v, which in our case is that of our 
LABOCA measurement, 345 GHz. 

The fit to the far-IR/sub-mm observational points gives us 
Td ~ 95 + 9 K andyS ~ 1.10 ± 0.14, which allows us to calculate 
a CSE dust mass of 0.010 + 0.003 Mq. The uncertainties listed 
are the statistical errors resulting from the fitting procedure and 
from error propagation in Eq.|2] 

^astro-Carrizo et al. (200lb derive a formula to estimate the 
mass of atomic gas from the flux observed in the [CII] line at 
157.8 jt/m, which we report here: 



M(Mo) = 7.0 X 10*^ F[c//](erg cm-^^"') D(kpc)2 FJXc 



(8) 



where F^cii] is the observed flux in the [CII] line, D the distance 
to the star, Xc the abundance of C ( 1.5x10"^ for C-rich stars and 
2x10 '* for O-rich sources), and Fc is a correction factor that 
accounts for low excitation temperatures. The correction factor 
is always > 1 and depends on the density and kinetic temperature 
of the gas, but at the densities found in proto-PNe it is unlikely 
to be larger than 2 or 3. 

From Eq. [8] we can then calculate the mass of atomic gas in 
the nebula around 115452 to be about 1.9 M© D(kpc)^, assuming 
Fc - I- 115452 appears to be surrounded by a massive CSE: 
still assuming Fc - 1, it would contain about 12 M© of gas, at 
the distance of 2.5 kpc. A statistical error on the mass estimate 
may be calculated from the error on the flux of the [CII] line, 
but this would not account for the real uncertainty on the mass, 
which is much larger and due to the degree of contamination. 
The high value of mass found is beyond what is expected from 
intermediate-mass stars and contrasts with the small momentum 



found above for the molecular component, which is considered 
as an indication for the source to be a low-mass star. Such a large 
estimate for the circumstellar atomic mass can be regarded as an 
indication for ISM contamination of the far-IR spectrum. 

For a gas-to-dust mass ratio of 200 - a standard value in 
evolved stars - our dust mass estimate implies a gas mass of 
about 2 Mo. Contamination would then enhance the flux in the 
[CII] line by a factor of 6. 



5. Discussion and conclusions 

We have observed a sample of stars where the detection of maser 
velocities larger than what observed on the AGB indicates that 
these sources have started their post-AGB evolution. In particu- 
lar, we have included four of the water-fountain nebulae known 
that were not previously investigated for CO emission. As our 
targets are located in the Galactic plane, spurious emission from 
the interstellar medium typically contaminates the spectra and 
does not allow us to obtain conclusive detections. By comparing 
the maser features with the CO lines, we find hints for possi- 
ble CO detections, but only with sensitive maps of the regions 
around the targets can a distinction between ISM and circum- 
stellar CO be made. 

Only IRAS 15452 - 5459 is clearly detected in both the 
J = 3 - 2 and J = 4 - 3 lines. The star shows in its central region a 
fast CO outflow with velocities of about Vexp - 45 km s"', while 
narrow peaks point to an expansion velocity of ~ 14 km s"'. 
The latter falls within the range of velocities observed during 
the AGB and, as it is detected at all pointings, it can be inter- 
preted as due to the remnant AGB wind. With its four peaks, 
the CO feature resembles what observed in the OH maser and 
points to a complex emitting structure, rather than an expand- 
ing spherical envelope. The broad component observed only in 
the central region of the source points instead to an outflow of 
molecular gas. If such an outflow is in the form of a collimated 
wind can be estimated only with interferometric observations. 
The mass-loss rate calculated for this component is quite high 
(~ 10""^ Mq yr"') but the assumed spherical symmetry implies 
that the value is likely overestimated. We have also calculated the 
scalar momentum and kinetic energy associated to the outflow- 
in g material, which we ha ve compared with the results obtained 
bv lBuiarrabal et all(l2001l) . concluding that radiation pressure on 
dust grains alone cannot explain the dynamics of the CSE. 
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Fig. 5. Spectra of the J = 3 - 2 line observed 10" North-East 
of the central source in IRAS 15452 - 5459 (top), on the cen- 
tral source (center) and 10" South- West of the central position 
(bottom). Overplotted in red are Gaussian fittings to the narrow 
features and a parabola for the broad component in the central 
pointing. 



By fitting a modified blackbody to the far-lR/sub-mm obser- 
vational points of 115452, we have estimated a dust mass in its 
CSE of ~0.01 M0, while through the analysis of its ISO spec- 
trum, we have calculated an envelope gas mass of about 12 M©. 
The uncertainty in the latter value is dominated by the degree of 
contamination of the [Cll] line. The derived values of gas and 
dust masses (12 and 2 Mq, respectively) indicate an enhance- 
ment factor of 6 of the line flux because of contamination, if we 
assume a standard value of 200 for the gas-to-dust mass ratio. 

Nevertheless, the total CSE mass of 115452 may be larger 
than 2 M© and a higher gas-to-dust mass ratio should be consi- 
de red. 115452 belon gs to a group of OH masing stars classified 
by ISevensterl (|2002|) as LI, because they lie on the left of the 
evolutionary sequence on the IRAS color-color diagram, while 
"standard" post-AGB stars lie on the right side (RI sources). LI 
stars have higher outflow velocities and mass-loss rates and they 
may loop back to the left of the evolutionary sequence because 




Fig. 7. Spitzer composite of the region around 115452 from the 
GLIMPSE and MIPSGAL archives. Images at 5.8, 8.0 and 24 
jum are coded on a logarithmic scale as blue, green and red, re- 
spectively. The image covers an area about 24' x 19' wide. North 
is up and East left. 
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Fig. 8. SED of 115452 with the ISO spectrum shown in red and 
the DUSTY model as a black solid line. The data from Table |6] 
and our measurement at 345 GHz (870 fj.m) are shown as solid 
circles. 



of an early AGB termination, which would also cau se them to 
host O H masers for a much longer time than R I stars (ISeyensten, 
l2002h . Finally, by scale height considerations. ISevensteJ (l2002h 
argue that LI stars may be the precursors of bipolar planetary 
nebulae and have larger masses than RI sources (ending in ellip- 
tical planetary nebulae), wit h average values of 4 and 1.7 Mq 
for LI and RI, respectively (iDeacon et al.L l2007h . Indeed, the 
bipolarity of the nebula around 115452 has been proven by the 
Hubble images. For this star to have a CSE mass within values 
adequate for LI sources (a minimum of 3-4 M©), its gas-to-dust 
mass ratio must be at least 300-400, on the basis of our dust 
mass estimate. This would decrease the degree of contamination 
of the [CII] line mentioned above to a factor of 3-4. 

115452 is certainly a peculiar object under many points of 
view. For example, the 4 peaks seen in its satellite OH line are al- 
most unique among the several hundreds of evolved stars where 
this feature has been detected. Another target with such a profile 
of its OH maser is OH 19.2 - 1.0, in which the spectral shape 
of the features shows a smaller degree of symmetry than that 
in 115452. The four peaks in 115452 are stunningly similar to 
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([01] 63 /xm + [CII] 158 /xm) / I„^ 
10° 



10^ 10^ 10^ 10^ 10^ 10^ 10^ 
n (cm"^) 

Fig. 9. Contour plots of the ratio of the intensity of the [CII] line 
plus that of the [OI] line over the FIR continuum as a function 
of the intensity of the FUV field G and the atomic H density n. 
The dashed contour corresponds to the value of 0.0037 found for 
115452 and the solid circle corresponds to G = 170 and n = 117. 
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one of the profiles derived theoretically by iGrinin & Grigor' evi 
(Il983h in the case of an unsaturated maser origina ting from 
a ro tating and expanding disk. As pointed out by Chapman 
(IT988) for OH 19.2-0.1, t he applicability of the calculations by 
iGrinin & Grigor'evl(ll983b to an evolved star is dubious, because 
OH masers in these sources are typically saturated. Nevertheless, 
the mid-IR excess of its SED and the mixed chemistry may ac- 
tually point to the existence of a disk in 115452. The detection of 
the PAH features in a region of the ISO spectrum that does not 
appear to be affected by ISM contamination and the presence of 
O-bearing molecules can in fact be explained by a circumstellar 
disk. Such disk would serve as a reservoir of O-rich material, 
where the masers would arise, while the PAHs would be located 
in the circumstellar envelope, converted into C-rich by the third 
dredge-up. 

Observations at high angular resolution of the mid-IR spec- 
trum, the masers and the CO outflow are necessary to reconstruct 
the morphology of the central emitting region in 115452. 
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